CALDERON'S INVERSE PROBLEM FOR ANISOTROPIC 
CONDUCTIVITY IN THE PLANE 



KARI ASTALA, MATTI LASSAS, AND LASSI PAlVARINTA 

Abstract: We study inverse conductivity problem for an anisotropic conductivity 
a G L°° in bounded and unbounded domains. Also, we give applications of the 
results in the case when Dirichlet-to-Neumann and Neumann-to-Dirichlet maps are 
given only on a part of the boundary. 

1. INTRODUCTION 
Let us consider the anisotropic conductivity equation in two dimensions 

(1) V . ffVu= j2^*(x)JLu = Oinfi, 

j,k=i 

u\dn = 4>- 

Here O C I 2 is a simply connected domain. The conductivity a = [c- ?fe ]| jfe=1 is a 
symmetric, positive definite matrix function, and <\> G H l l 2 (dVl) is the prescribed 
voltage on the boundary. Then it is well known that equation |p has a unique 
solution u G H 1 ^). 

In the case when a and dfl are smooth, we can define the voltage-to-current (or 
Dirichlet-to-Neumann) map by 

(2) K{<j>) = Bu\dn 
where 

(3) Bu = v- crVu, 

u G H 1 ^) is the solution of (jlj, and v is the unit normal vector of dfl. Applying 
the divergence theorem, we have 

(4) Q*M) ■= I E ^ x )^^ dx = J m A M<i>ds } 

where dS denotes the arc lenght on dQ. The quantity Qo-,q(0) represents the power 
needed to maintain the potential on <9f2. By symmetry of A CT , knowing Q a ,n is 
equivalent with knowing A a . For general Q and a G L°°(Q), the trace u\qq_ is defined 
as the equivalence class of u in H l {yt)/Hl{yi) (see 0) and formula JH) is used to 
define the map A CT . 
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If F : Q — > Q, F(x) = (F 1 (x), F 2 (x)), is a diffeomorphism with F\gn = Identity, 
then by making the change of variables y = F(x) and setting v = u o F^ 1 in the 
first integral in (JH), we obtain 



V- [F*<r)Vv = in ft, 



where 



2 OF* , , dF k 



det[fg-(x)] ^daiP K 1 dx^ J K J 



x=F~ 1 (y) 



or 



(6) = -L-DF{x) a(x) DF(x] 

is the push-forward of the conductivity a by F. Moreover, since F is identity at dQ, 
we obtain from that 

A Ft(T = A CT . 

Thus, the change of coordinates shows that there is a large class of conductivities 
which give rise to the same electrical measurements at the boundary. 

We consider here the converse question, that if we have two conductivities which 
have the same Dirichlet-to-Neumann map, is it the case that each of them can be 
obtained by pushing forward the other. 

In applied terms, this inverse problem to determine a (or its properties) from 
A CT is also known as Electrical Impedance Tomography. It has been proposed as a 
valuable diagnostic, see [TT] . 

In the case where a^ k (x) = cr(x)5 jfc , a(x) G M+, the metric is said to be isotropic. 
In 1980 it was proposed by A. Calderon p] that in the isotropic case any bounded 
conductivity cr(x) might be determined solely from the boundary measurements, i.e., 
from A CT . Recently this has been confirmed in the two dimensional case (c.f. [6j). In 
the case when isotropic a is smoother than just a L°°-function, the same conclusion 
is known to hold also in higher dimensions. 

The first global uniqueness result was obtained for a C°°-smooth conductivity in 
dimension n > 3 by J. Sylvester and G. Uhlmann in 1987 [36j. In dimension two 
A. Nachman [SH] produced in 1995 a uniqueness result for conductivities with two 
derivatives. The corresponding algorithm has been successfully implemented and 
proven to work efficiently even with real data [33} I27|. The reduction of regularity 
assumptions has since been under active study. In dimension two the optimal ir- 
regularity was obtained in [6j. In dimension n > 3 the uniqueness has presently 
been shown for isotropic conductivities a G W /3//2 '°°(f2) in jHI] and for globally 
C 1+£ -smooth isotropic conductivities having only co-normal singularities in [T3] . 

Also, the stability of reconstructions of the inverse conductivity problem have been 
extensively studied. For these results, see [21 13 IH where stability results are based 
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on reconstruction techniques of [8] in dimension two and those of [28] in dimensions 
n > 3. 

In anisotropic case, where a is a matrix function and the problem is to recover the 
conductivity a up to the action of a class of diffeomorphisms, much less is known. In 
dimensions n > 3 it is generally known only that piecewise analytic conductivities 
can be constructed (see |THl ) - For Riemannian manifolds this kind of technique 
has been generalized in |22[ I2*U1 121) . In dimension n = 2 the inverse problem has 
been considered by J. Sylvester [HE] for C 3 and Z. Sun and G. Uhlmann [31] for 
W /1 ' p -conductivities. The idea of [3E] and [31] is that under quasiconformal change 
of coordinates (cf. [HE]) an y anisotropic conductivity can be changed to isotropic 
one, see also section OH below. The purpose of this paper is to carry this technique 
over to the L°°-smooth case and then use the result of [H] to obtain uniqueness up 
to the group of diffeomorphisms. 

The advantage of the reduction of the smoothness assumptions up to L°° does 
not lie solely on the fact that many conductivities have jump-type singularities 
but it also allows us to consider much more complicated singular structures such as 
porous rocks [10]. Moreover it is important that this approach enables us to consider 
general diffeomorphisms. Thus anisotropic inverse problems in half-space or exterior 
domains can be solved simultaneously. This will be considered in Section [2] 

If O C I 2 is a bounded domain, it is convenient to consider the class of matrix 
functions cr = [a^ h ] such that 

(7) [a ij ] G L°°(ft;M 2x2 ), [a**]' = [a% C^I < [a ij ] < C I 

where C > 0. In sequel, the minimal possible value of C is denoted by C (a). We 
use the notation 

= {a G L°°{Q; R 2x2 ) | C (a) < oo}. 

Note that it is necessary to require C (o") < oo as otherwise there would be coun- 
terexamples showing that even the equivalence class of the conductivity can not be 
recovered [TH UK] . 

Our main goal in this paper is to show that an anisotropic L°°-conductivity can 
be determined up to a W /1,2 -diffeomorphism: 

Theorem 1. Let Q C M 2 be a simply connected bounded domain and a G L°°(Q; R 2x2 ). 
Suppose that the assumptions (GJ) are valid. Then the Dirichlet-to- Neumann map A CT 
determines the equivalence class 

E a = {cr 1 G E(O) | o~i = F*a, F : Q — ► Q is W 1,2 -diffeomorphism and 

F\dn = I}. 

We prove this result in Sectional 

Finally, note that the iy 1,2 -diffeomorphisms F preserving the class E(fi) are pre- 
cisely the quasiconformal mappings. Namely, if cr G E(O) and U\ = F*(a ) G 
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E(F(Q)) then 

(8) i|DF(a;)|| 2 / < DF(x) a (x) DF{x) t < dJ F (i)/ 
Co 

where I = and we obtain 

(9) \\DF(x)\\ 2 < KJ F (x), fora.e. xett 

where K = CiC < oo. Conversely, if © holds and F is PF^' 2 -homeomorphism then 
F*a G E(F(f2)) whenever o G E(fi). Furhtermore, recall that a map F : Q — > Q 
is quasiregular if F G W l( £ (O) and the condition (jUJ) holds. Moreover, a map F is 
quasiconformal if it is quasiregular and a W /1 ' 2 -homeomorphism. 

2. CONSEQUENCES AND APPLICATIONS OF THEOREM □ 

Here we consider applications of the diffeomorphism-technique to various inverse 
problem. The formulated results, Theorems are proven in Section HJ 

2.1. Inverse Problem in the Half Space. Inverse problem in half space is of cru- 
cial importance in geophysical prospecting, seismological imaging, non-destructive 
testing etc. For instance, the imaging of soil was the original motivation of Calderon's 
seminal paper [H]. As we can use a diffeomorphism to map the open half space to 
the unit disc, we can apply the previous result for the half space case. One should 
observe that in this deformation even infinitely smooth conductivities can become 
non-smooth at the boundary (e.g. conductivity oscillating near infinity produces a 
non-Lipschitz conductivity in push-forward) and thus the low-regularity result [H] is 
essential for the problem. 

Thus, for a G E(R 2 ) let us consider the problem 

(10) V-crVu = in R 2 _ = {(x\x 2 ) \x 2 < 0}, 

(11) u\ dR 2_ = 0, 

(12) u G L^iRl). 

Notice that here the radiation condition at infinity (fl2|l is quite simple. We assume 
just that the potential u does not blow up at infinity. The equation (ffnTTl2|) is 
uniquely solvable and as before we can define 

K : HUl p {dR 2 _ ) - H'^idRi), <j>^v aVu\ ml . 

Theorem 2.1. The map A CT determines the equivalence class 

E a = {tri G E(IR 2 ) | a 1 = F*a, F : R 2 _ -> R 2 _ is W 1 ' 2 -diffeomorphism, 

F\dRl — 

Moreover, each orbit E a contains at most one isotropic conductivity, and conse- 
quently if a is known to be isotropic, it is determined uniquely by A a . 
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Note that the natural growth requirement lim| 2 |_ >00 = oo follows automati- 

cally from the above assumptions on F. 

2.2. Inverse Problem in the Exterior Domain. An inverse problem similar to 
that of the half space can be considered in an exterior domain where one wants to 
find the conductivity in a complement of a bounded simply connected domain. This 
type of problem is encountered in cases where measurement devices are embedded 
to an unknown domain. 

In the case of S = M. 2 \ D, where D is a bounded Jordan domain, we consider the 
problem 

(13) V- aVu = in S, 

(14) u\ 9s = 0Gtf 1/2 (aS), 

(15) u G L°°(S). 

Again, the radiation condition (fT5j) of infinity is only that the solution is uniformly 
bounded. For this equation we define 

K : H^idS) -> H- l l\dS), <f>^u- aX7u\ as . 

Surprisingly, the result is different from the half-space case. The reason for this is 
the phenomenon that the group of diffeomorphisms preserving the data does not 
fix the point of the infinity. More precisely, there are two points x ,x\ G S U {oo} 
such that F(xq) = oo, F _1 (a; 1 ) = oo, and F : S \ {x } — > S \ {xi}. In particular, 
this means that the uniqueness does not hold up to diffeomorphisms mapping the 
exterior domain to itself. 

For convenience, we compactify S by adding one infinity point, denote S — S U 
{oo}, and define cr(oo) = 1. We say that F : S — > S is a W /1 ' 2 -diffeomorphism if F 
is homeomorphism and a W /1 ' 2 -diffeomorphism in spherical metric 

Theorem 2.2. Let a G S(S'). Then the map A a determines the equivalence class 

E a ,s — Wi G I 0i = F*a, F : S —> S is a W 1 ' 2 -diffeomorphism, 

F\ dS = I}. 

Moreover, if a is known to be isotropic, it is determined uniquely by A CT . 

2.3. Data on Part of the Boundary. In many inverse problems data is measured 
only on a part of the boundary. For the conductivity equation in dimensions n > 3 
it has been shown that if the measurements are done on a part of the boundary, then 
the integrals of the unknown conductivity over certain 2-planes can be determined 
[T2] . In one-dimensional inverse problems partial data is often considered with two 
different boundary conditions, see e.g. j24j[2H]. For instance, in the inverse spectral 
problem for a one-dimensional Schrodinger operator, it is known that measuring 
spectra corresponding to two different boundary conditions determine the potential 
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uniquely. Here we consider similar results for the 2-dimensional conductivity equa- 
tion assuming that we know measurements on part of the boundary for two different 
boundary conditions. 

Let us consider the conductivity equation with the Dirichlet boundary condition 

(16) V-aVu = Oinfi, 

u\dn = 4> 
and with the Neumann boundary condition 

(17) V-aVv = OmQ, 

v-aVv\m = ip, 

normalized by J dn vdS = 0. Let T C dfl be open. We denote by Hq(T) the 
space of functions / G H s (dQ) that are supported on T and by H S (T) the space of 
restrictions /|r of / G H s (dQ). We define the Dirichlet-to-Neumann map Ar and 
Neumann-to-Dirichlet map Er by 

A r : tf 1/2 (r) - H-^iT), ^ {v aVu)\ T , 
E r : H- 1/2 (T) - H X /\T), V » «lr- 

Theorem 2.3. Let T C dQ be open. Then knowing dQ and both of the maps Ar 
and Er determine the equivalence class 

E a ,r = ^ E(f2) | ox = F^a, F : Q ^ Q is a W 1,2 -diffeomorphism, 

F\ r = I}. 

Moreover, if a is known to be isotropic, it is determined uniquely by A-p and Er. 

3. PROOF OF THEOREM □ 

3.1. Preliminary Considerations. In the following we identify M 2 and C by the 

map (x 1 , x 2 ) I— > x 1 + ix 2 and denote z = x 1 + ix 2 . We use the standard notations 

d z = ~(<9i -id 2 ), dz = ^(d 1 +id 2 ), 

where dj = djdxK Below we consider o : Q — > M 2x2 to be extended as a function 
o \ C — > M 2x2 by defining a(z) = I for 2; G C\f2. In following, we denote Co = Cq(ct). 
For the conductivity a = a^ k we define the corresponding Beltrami coefficient (see 

geheied 

-q^ + a 22 ^)-^ 12 ^) 

(18) Mil- 2 ) = ; -■ 

' a u {z) + a 22 {z)+2^/det{cx{z)) 

The coefficient satisfies |/ii(z)| < k < 1 and is compactly supported. 

Next we introduce a W /1 ' 2 -diffeomorphism (not necessarily preserving the bound- 
ary) that transforms the conductivity to an isotropic one. 
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Lemma 3.1. There is a quasiconformal homeomorphism F : C — > C such that 

(19) F(z) = z + 0(^) as \z\ -> oo 

and such that F e W^(C; C),2<p< p(C ) = for which 

(20) {F*a){z) = a{z) : = det(a(F^(z)))^ . 

Proof. The proof can be found from for C 3 -smooth conductivities, see also 
[T7] . Because of varying sign conventions, we sketch here the proof for readers 
convenience. We need to find a quasiconformal map F such that 



(21) DF a DF 1 = y/det(a) J F I 

where Jf = det(DF) is the Jacobian of F. Denoting by G = [<7y]fj = i the inverse of 
the matrix a / \J det(cr) we see that the claim is equivalent to proving the following: 
For any symmetric matrix G with det(G) = 1 and < G < Kl there exists a 
quasiconformal map F such that 

(22) J F G = DF l DF. 

Next, the non-linear equation (|22|) can be replaced in complex notation by a linear 
one. Indeed, if F = u + iv then (|22|) is equivalent to 



(23) Vv = JG -1 Vu, where J = ^ J ^ ^ 

This follows readily from the identity 

DF 1 ,] = det(DF) J (DF)- 1 = JG^DF 1 



where the latter equality uses (J22J). The matrix J corresponds to the multiplication 
with the imaginary unit % in complex notation. Denoting by R = ( J (the 



v -1 f 

matrix corresponding to complex conjugation) we see that (|23|) is equivalent to 

(24) Vu + JVv = (G - 1)(G + ly^Vu - JVv). 

But, Vu + JVv = 2c\F and R(Vu — JVv) = 2d z F in complex notation and hence 
(|2~H) becomes 

(25) d?F = tn(z)d z F 
where 

fjL 1 = (G- 1)(G + l)" 1 ^ = (Vdet a I - a)(VdetaI + a^R. 
A direct calculation gives 
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which shows that the matrix [i\ = (G — 1)(G + 1)~ 1 R corresponds to a multiplication 
operator (in complex notation) by the function 

/ \ _ 9u {*) ~ 922(z) + 2ig 12 (z) 
2 + ^11(2) +6-22(2;) 

This gives (|T£|) since G^ 1 = a / a/ det(a). Since \fii(z)\ < k < 1 for every z G C it is 
well known by [TJ Thm. V.l, V.2] that the equation (OKI) with asymptotics 

F(z) = z + 0(~), as z -> 00 

z 

has a unique (quasiconformal) solution F. The fact that F G W lt £ (C; C), 2 < p < 
^ follows from [5J. □ 

In this section we denote by F = F a the diffeomorphism determined by Lemma 
13.11 We also denote fl = F(Q) where F is as in Lemma l3~Tl Note that (fT9|) implies 
also that 

(26) F-\z) = z + 0(^) as \z\ ^00. 

Later we will use the obvious fact that the knowledge of map A CT is equivalent to 
the knowledge of the Cauchy data pairs 

C a = {Han, v oVu\ m ) \ u G H\Q), V- aVu = 0}. 

In addition to the anisotropic conductivity equation JTJ we consider the corre- 
sponding conductivity equation with isotropic conductivity. For these considera- 
tions, we observe that if u satisfies equation (JTJ) and a is as in (|2T)|) then the function 

w(x) =u(F- 1 (x)) G H\h) 

satisfies the isotropic conductivity equation 

(27) V-?Vw = infi, 

w \dh = <P°F-\ 

Thus, a can be considered as a scalar, isotropic L°°-smooth conductivity aL We 
continue also the function a : Q — > M. + to a function a : C — > M. + by defining 
a(x) = 1 for x G C \ Q. 

3.2. Conjugate Functions. While solving the isotropic inverse problem in [6J, the 
interplay of the scalar conductivities a(x) and played a crucial role. Motivated 
by this, we define 

det{a{x)) 

Note that for a isotropic conductivity a = 1/er. 



CALDERON'S INVERSE PROBLEM FOR ANISOTROPIC CONDUCTIVITY IN THE PLANE 9 

Let now F be the quasiconformal map defined in Lemma I3~T1 and a = F*a as in 
(f20|) . We say that w G i? 1 (SI) is a a-harmonic conjugate of w if 

(28) d!w(z) = -a(z)d 2 w(z), 

d2~w(z) = a{z)d\w{z) 

for z = x 1 + ix 2 G C. Using w we define the function u that we call the a-harmonic 
conjugate of u, 

u(x) = w(F(x)). 
To find the equation governing u, it easily follows that (c.f. [£]]) 

(29) V- hvw = in O, 

a 

and by changing coordinates to y = F(x) we see that l/er = Fjj. These facts imply 

(30) V- aVu = in Q. 

Thus u is the a-harmonic conjugate function of u and we have 

(31) Vw = JaVu, Vw = JaS/u. 

Since u is a solution of the conductivity equation if and only if u + c, c G C, is 
solution, we see from (|3TJ) that the Cauchy data pairs C a determine the pairs C% 
and vice versa. Thus we get, almost free, that A CT determines A^, too. 
Let us next consider the function 

(32) f(z)=w(z)+iw(z). 

By 0, it satisfies the pseudo-analytic equation of second type, 



(33) dgf = Hz d z f 
where 

/ \ ~ , l — a(z) ,~ , si Co — 1 

(34) ^ 2 { Z ) = ^J- ^(z < < i. 

l + a(z) G + 1 

Using this Beltrami coefficient, we define // 2 = ° F. 
We will need the following: 

Lemma 3.2. Let g = f o F where F : ^ is a quasiconformal homeomorphism 
and f is a quasiregular map satisfying 



(35) dzf = V-idzf and d^F = md z F, 

where JI2 = 1^2 F and \i\ satisfies \fij\ < K < 1 and j2 2 is real. Then g is 
quasiregular and satisfies 

(36) ckg = v x d z g + v 2 d z g, 
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where 

1 — /in 1 — |/il| 2 

(37) v\ = /ii- ; — 7^-5, and v 2 = /U 2 - , — tz-$. 

Conversely, if g satisfies pity) with v 2 real and \v\\ + \u 2 \ < k! < 1 then there exists 
unique (i\ and // 2 such that (TF/p holds and f — g o F" 1 satisfies pity) . 

Proof. We apply the chain rule 

d(foF) = (df)oF-dF+(df)oF-dF, 
d(foF) = (df)oF-dF+(df)oF-dF, 

and obtain 



v x d z g + v 2 d z g — dfoF- dF- {v x + z^Z-ti/^) + df o F- dF- (z/ 2 + vx^ 2 ) 

and 



dgg = ny df o F- dF + /i 2 - 9/ o F- dF. 
Hence, if /ii,/i 2 , z/ 1? and z/ 2 are related so that 

fJL X = V X + V1H2, tl2 = V 2 + F x /i2, 

we see that (f36|) and (|3TJ) are satisfied. 

It is not difficult to see that for each V\ and v 2 (J3TJ) has a unique solution \i\ y \i 2 
with < «' < 1, j — 1,2. Again, the general theory of quasiregular maps [l] 
implies that (|35|) has a solution and the factorization g = f o F holds. □ 

Note that lj3*7J) implies that 
38 + z/ 2 = TT\ — il — i - TT~^ < L 

I + I/X1M/X2I 1 + K 

Lemma 13.21 has the following important corollary, that is the main goal of this 
subsection. 

Corollary 3.3. If u £ H 1 ^) is a real solution of the conductivity equation (QJ), 
there exists u G H 1 ^), unique up to a constant, such that g = u + iu satisfies $36}) 
where 

, QO , o*-<P-1vP . l-det(a) 

(39) V\ = — t— r-, and v 2 



1 + tra + det(a) 1 + tr a + det(cr) 

Conversely, if V\ and v 2 , \v\\ + \v 2 \ < k' < 1 are as in Lemma \S.S\ then there are 
unique a and a such that for any solution g of $36\) u = Reg and u = Img satisfy 
the conductivity equations 

(40) V- aVu = 0, and V- aWu = 0. 
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Proof. Since g = foF where f = w + iw according to (jH2J), we obtain immediately 
the existence of u = w o F. Thus we need only to calculate i>\ and v 2 in terms of a. 
Note that by (JTHJ) , 

- tr(a)-2det((7) 1 / 2 



(41) M 2 = 
We recall that 

(42) fi2 
and thus 



tr(cr) + 2det(a) 1 /2' 
l-detfa) 1 ^ 



1 - \fii\ 2 f4 



1 + det(o-) 1/2 
4(det(a) 1 / 2 tr (a) + (1 + det(a)) det(a) 1 / 2 ) 



(1 + det(a) 1 /2)2( tr ( a ) + 2 det(a) 1 /2) 

which readily yields (|39|) from (j3"TJl . 

Note that since v\ and v 2 uniquely determine \i\ and /x 2 , they by (J4*T| and (f4*2|) 
also determine det(a) and tr(cr). After observing this, it is clear from (|39|) that o is 
uniquely determined by V\ and z/ 2 . □ 

Now one can write equations (|31|) in more explicit form 

(43) r- Vu\ a u = A a (u\ aQ ) 

where r = (—v 2 , v\) is a unit tangent vector of dfl. As Re g\gn = u\qq and Imgl^ = 
u\dQ,i we see that A CT determines the cx-Hilbert transform Ti a defined by 

(44) n a : H^ 2 (dn) -> H^idty/c, 

Put yet in another terms, for u,u E H l l 2 (dVL), u = 7i a u if and only if the map 
9(C) — ( u + £ e extends to f2 so that ((BE} is satisfied. 

Summarizing the previous results, we have 

Lemma 3.4. The Dirichlet-to- Neumann map A CT determines the maps and Ha- 

3.3. Solutions of Complex Geometrical Optics. Next we consider exponen- 
tially growing solutions, i.e., solutions of complex geometrical optics originated by 
Calderon for linearized inverse problems and by Sylvester and Uhlmann for non- 
linear inverse problems. In our case, we seek solutions G(z,k), z G C \ fl, i 6 C 
satisfying 

(45) dgG(z, k) = for z G C \ H, 

(46) G(z,k) = e lkz (l + O k (^)), 

(47) lmG(z,k)\ z£dn =n a (ReG(z,k)\ zedn ). 

Here Ok(h(z)) means a function of (z, k) that satisfies \Ok{h{z))\ < C(k)\h(z)\ for 
all z with some constant C{k) depending only on k G C. For the conductivity a we 
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consider the corresponding exponentially growing solutions W(z,k), z G C, k G C 
where 



(48) dgW(z, k) = Jl 2 (z)d z W(z, k), for z G C, 

(49) W(z,k) = e**(l + O k (±)). 

Note that in this stage, z G(z, k) is defined only in the exterior domain C\Q but 
z i— > W(z, k) in the whole complex plane. These two solutions are closely related: 

Lemma 3.5. For all k G C we have: 

i. The system HJffy has a unique solution W(z, k) in C. 

ii. The system f^ffT^ has a unique solution G(z, k) in C \ Q. 

Hi. For z G C \ Q we have 

(50) G(z,k) = W(F(z),k). 

Proof. For the claim i. we refer to [HI Theorem 4.2]. 

Next we consider ii. and iii. simultaneously. Assume that G(z, k) is a solution of 
(j43H4Tjl . By Lemma E21 and boundary condition (|4Tj) we see that equation 



(51) H{z,k) 



dh(z, k) = V\d z h + is 2 d z h } in fi, 

h(-,k)\ dn = G(-,k)\ d n 

has a unique solution where z/i and v 2 are given in (JHEJ). 
Let 

'G{z,k) fovzeC\n 
h(z, k) for zGO 

and H(z, k) = H(F~ 1 (z), k). Then H(z, k) satisfies equations 

dgH(z, k) = 0, for zGC \ fi, 

dzH(z,k) = jj, 2 (z)d z H(z,k), for zGO, 

and traces from both sides of dfl coincide. Thus H(z, k) satisfies equation in (jIBJ). 
Now (JH and flU yield that 

(52) H(z,k) = H{F~ 1 {z), k) 

= e X p(ikF- 1 (z))(l + O k ( 1+ ^ )) 

= exp(iA;z)(l + e> ifc (— -- )) 

showing that if satisfies (pl8Tli9"|) . Thus by i., H(z,k) = W(z, k). This proves both 
ii. and iii. □ 
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3.4. Proof of Theorem [H As G(z, k) is the unique solution of (I3HH37}) and the 
operator appearing in boundary condition (fTTj) is known, Lemma,ta l3.4l a,nri l3.5l impry 
the following: 

Lemma 3.6. The Dirichlet-to-Neumann map A a determines G(z,k), z £ C \ Q, 
keC. 

Next we use this to find the diffeomorphism F a outside Q. 

Lemma 3.7. The Dirichlet-to-Neumann map A CT determines the values the restric- 
tion F a \c\n- 

Proof. By (jHOjl . G(z, k) = W(F(z), k), where W(z, k) is the exponentially growing 
solution corresponding to the isotropic conductivity a. Thus by applying the sub- 
exponential growth results for such solutions, [6, Lemma 7.1 and Thm. 7.2], we have 
representation 

(53) W(z, k) = exp(ikip(z, k)) 
where 

(54) lim sup \<f(z, k) — z\ = 0. 

As F(z) = z + D(l/z), and G(z, k) = W{F(z), k) we have 

(55) lim lo S G ( z ' fc ) = lim v (F(z), k) = F(z). 

k— >oo tk k^oo 

By Lemma l3~6l we know the values of limit (|55| for any z £ C \ Q. Thus the claim 
is proven. □ 

We are ready to prove Theorem [TJ 
Proof. As we know F|c\n £ W 1,p , 2 < p < p(C ), we in particular know Q = 
C \ (F(C \ ft)). When u is the solution of conductivity equation JTJ with Dirichlet 
boundary value and w is the solution of (|2*7jl with Dirichlet boundary value = 
4> o h, where h = F -1 ]^ we see that 



(56) / 4>^<t>dS = Q~ h {w) = Q^{u)= <j>K4>dS. 

Jan Jan 

Here, the second identity is justified by the fact that F is quasiconformal and hence 
satisfies ©. Since A CT and A^ are symmetric, this implies 



(57) / if>Az<pdS= / ipA^dS 

Jan Jan 

for any ip,(f) £ ^^(dVL) and ip,(f) £ H l l 2 (dVt) are related by ip = ipoh and ip = ipoh. 
Note that £ H l ! 2 (dti) if and only if = o h £ H l / 2 {dVt). To see this, extend 
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to a if function and after that define <fi in the interior of fl by (ft = (p o F 1 . 
Now 

1 1 W| |i» ( n) ~ J V0- aVjdx~ Vf ^4>dx ~ 1 1 V0| |* 2{S) 
and hence 

ll^llHVa ( an) ~ ll^llffi(n) ~ H^l^) ~ H^l/n/ a (an)- 

As we know F|c\n and A CT , we can find A^ using formula ijHTJ) . By 0, the map 
Asp determines uniquely the conductivity a on ft in a constructive manner. 

Knowing ft, ft, and the boundary value / = F\qq of the map F : ft — > f2, we next 
construct a sufficiently smooth diffeomorphism if : f2 — ► fi. First, by the Riemann 
mapping theorem we can map fl and fl to the unit disc D by the conformal maps 
R and R, respectively. Now 

G = Ro F o R- 1 : D ^ D 

is a quasiconformal map and since we know R and i?, we know the function g = G\gn 
mapping <9D onto itself. The map g is quasisymmetric (cf. [T]) and by Ahlfors- 
Beurling extension theorem (TJ Thm. IV. 2] it has a quasiconformal extension AB(g) 
mapping D onto itself. Note that one can obtain AB(g) from g constructively by 
an explicit formula |Tj. p. 69]. Thus we may find a quasiconformal diffeomorphism 
H = R- 1 o lAB(g)}- 1 oR,H:fl^fl satisfying H\ dh = 

Combining the above results, we can find HJj that is a representative of the 
equivalence class E a . □ 

In the above proof the Riemann mappings can not be found as explicitly as the 
Ahlfors-Beurling extension. However, there are numerical packages for approxima- 
tive construction of Riemann mappings, see e.g. |26j. 



4. PROOFS OF CONSEQUENCES OF MAIN RESULT 
Here we give proofs of Theorems I2.1H2.31 

Proof of Theorem EjJ Let F : R 2 _ = R + zIR_ -> D be the Mobius transform 



Since this map is conformal, we see that C (-F*cr) = C (a). Let a = F*a be the 
conductivity in D. Then K^<p is determined as in l(57jl for all <fi G C~(<9IB) \ {1}). 

Since A 5 1 = and functions C © C^(dB \ {1}) are dense in the space H 1 / 2 (dB), 
we see that A CT determines the Dirichlet-to-Neumann map A^ on dU>. Thus we can 
find the equivalence class of the conductivity on O. Pushing these conductivities 
forward with F~ l to R?., we obtain the claim. □ 
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Proof of Theorem 12. 2L Let F : S — > O \ {0} be the conformal map such that 

lim F(z) = 0. 

z— >oo 

Again, since this map is conformal we have for a = F*a the equality Cq(ct) = 
Co(F*a). Moreover, if u is a solution of (fT3j) . we have that to = aoF -1 is solution of 

(58) V-?Vw = 0inO\{0}, 

w\do = 4 > °F~ 1 , 
w e L°°(D). 

Since set {0} has capacitance zero in D, we see that u> = W|o\{o} where 

(59) V- ctVW 7 = in D, 

W|eD = (j)oF'\ 

Since F can be constructed via the Riemann mapping theorem, we see that A CT 
determines A^ on <9D and thus the equivalence class E„. When F : D — > D is a 
boundary preserving diffeomorphism, we see that F _1 o F o F defines a diffeomor- 
phism S — > S. Since we have determined the conductivity a up to a boundary 
preserving diffeomorphism, the claim follows easily. □ 

Proof of Theorem [231 Let © C C be the unit disc and D + = {z G © | Re z > 

0}. Let F : Q — > D + be a Riemann mapping such that 

D + cRxl + , F(T) = <9D + \(Rx {0}), F(dtt \ T) = <9B+ n (R x {0}). 
Let rj : (a; 1 ,^ 2 ) t— > (x 1 , —x 2 ) and define D_ = ?y(D + ), and a = F^a. Let 




for x G D+, 
for x G D_. 



Consider equation 

(60) V-dVw = OinD. 

Using formula (|5TJ) we see that F and A r determine the corresponding map Ap(r) 
for a. Similarly, we can find ^f(v) f° r 

Then A_p(r) determines the Cauchy data on the boundary for the solutions of (jfiUJ) 
for which w G if 1 (ED), w = —wot]. On the other hand, £_F(r) determines the Cauchy 
data on the boundary of the solutions of (|60|) for which w G if 1 (D) and w = w o rj. 
Now each solution w of (f60|) can be written as a linear combination 

1 1 

w(x) = ~{w(x) + w(t}(x))) + —(w(x) — w(r](x))). 

Thus the maps A_F(r) and £,F(r) together determine C$, and hence we can find a up 
to a diffeomorphism. We can choose a representative <Jo of the equivalence class Fj 
such that <7o = <?o ° V- I n f ac t, choosing a symmetric Ahlfors-Beurling extension in 
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the construction given in the proof of Theorem HJ we obtain such a conductivity. 
Pushing the conductivity &o from D + to Q with we obtain the claim. □ 
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